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miR-16 and miR-26a target checkpoint kinases Wee1 
and Chk1 in response to p53 activation by genotoxic 
stress 

L Lezina^'^'^, N Purmessur^'^, AV Antonov^, T Ivanova^, E Karpova^, K Krishan^, M Ivan^, V Aksenova\ D Tentler\ AV Garabadgiu\ 
G Melino^'^'^ and NA Barlev*'^'^'^ 

The tumour suppressor p53 is a crucial regulator of cell cycle arrest and apoptosis by acting as a transcription factor to regulate 
a variety of genes. At least in part, this control is exerted by p53 via regulating expression of numerous microRNAs. We identified 
two abundantly expressed microRNAs, miR-16 and miR-26a, whose expression is regulated by p53 during the checkpoint arrest 
induced by the genotoxic drug, doxorubicin. Importantly, among the targets of these mlRs are two critical checkpoint kinases, 
Chk1 and Wee1. The p53-dependent augmentation of miR-16 and miR-26a expression levels led to the cell cycle arrest of tumour 
cells in G1/S and increased apoptosis. Strikingly, the bioinformatics analysis of survival times for patients with breast and 
prostate cancers has revealed that co-expression of mir-16 and miR-26a correlated with a better survival outcome. Collectively, 
our data provide a novel mechanism whereby p53 represses Chkl and Weel expression, at least partially, via upregulation of 
miR-16 and miR-26a and thus sensitizes tumour cells to genotoxic therapies. 
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TP53 is a universal tumour suppressor whose major function 
is to control the stability of genome in response to various 
forms of genotoxic stress. Acting mostly as a transcription 
factor, p53 regulates the gene expression programme to 
prevent damaged cells from propagation. Products of 
p53-dependent genes induce DNA repair, cell cycle arrest 
in G1/S and G2/M phases, and apoptosis.'' G1/S arrest is 
mediated mostly by the p21 protein whose expression is 
induced by pSS.^""^ The former binds to and inhibits the 
kinase activity of cdk2/cyclin E and cdk4/cyclin D complexes, 
which are essential for G1/S transition. The p53-mediated 
G2 arrest involves, in addition to the p21 protein, other p53 
targets such as members of the 14-3-3 family of proteins^'^ 
and GADD45.'^ The 14-3-3 protein binds to and inhibits 
the phosphatase activity of Gdc25, whereas GADD45 
attenuates the activity of mitotic cdc2 (cdk1) kinase by 
sequestering its cyclin partner, cyclin B.^ Activity of cdkl is 
also regulated by two checkpoint kinases Ghekl and Weel . 
On DNA damage, Chkl and Weel phosphorylate and inhibit 
cdk1, thus preventing entry into mitosis.^"'''' CHK1 also 
phosphorylates p53 and together with ATM/ATR stabilizes 
and activates p53, further enhancing the G2-M checkpoint in 
the cell cycle. ""^ 

Recently, it has been demonstrated that microRNAs 
(miRNAs) are the integral part of the p53 network.^^'^^ 



miRNAs are a group of small (20-25 nucleotides) non- 
coding RNA molecules, which are expressed endogenously 
in cells and that regulate gene expression at the post- 
transcriptional level. To date, there is more than a 
dozen of different families of microRNA genes controlled 
by p53.''^'''^ Importantly, p53 has been shown to regulate 
both transcription and maturation of miRNA genes, whose 
products aid in the processes of cell cycle arrest and 
apoptosis. 

Here, we report an identification of two abundantly 
expressed microRNAs, miR-16 and miR-26a, whose 
expression is regulated by p53 during the genotoxic 
stress. Importantly, among the targets of these miRs are 
two critical checkpoint kinases, Chkl and Weel. 
The p53-dependent augmentation of miR-16 and miR-26a 
expression levels led to the cell cycle arrest of tumour 
cells in G1/S and increased apoptosis. Strikingly, the 
bioinformatics analysis of survival times for patients 
with breast and prostate cancers has revealed that 
co-expression of miR-16 and miR-26a correlated with a 
better survival outcome. Collectively, our data provide a 
novel mechanism whereby p53 represses Chkl and Weel 
expression, at least in part, via upregulation of miR-16 and 
miR-26a and thus sensitizes tumour cells to genotoxic 
therapies. 
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Results 

p53 and genotoxic stress control the expression of 
miR-16 and miR-26a. To define the spectrum of microRNA 
regulated by p53 on genotoxic stress, we employed the 
microarray expression analysis of microRNAs in a matching 
pair of cell lines that differ in their p53 expression 
status. Control U2-0S cells stably expressing non-specific 
scrambled shRNA (U2-0S scr) showed low levels of p53, as 
these cells overexpress p53-specific E3 ubiquitin ligase, 
Mdm2, which targets p53 for the 26S proteasome-dependent 
degradation. ""^"^^ A matching cell line (U2-0S KDp53) with 
suppressed expression of p53 due to the presence of 
p53-specific shRNA also showed very low level of p53 
expression. However, doxorubicin treatment stabilized p53 in 
U2-0S scr cells (Figure 1, left), whereas in U2-0S KDp53 
cells the level of p53 was still low (Figure 1a). 

To further assess the validity of the cell system used in our 
experiments, both U2-0S scr and U2-0S KDp53 cell lines 
treated with doxorubicin were analysed for their ability to elicit 
G 1 /S checkpoint arrest characteristic of the intact p53 function 
(Figure 1b). On genotoxic stress, U2-0S scr (p53 + ) cells 
exhibited cell cycle arrest in G1/S phase along with the 
accumulation of cells in G2/M. On the contrary, doxorubicin- 
treated U2-0S KDp53 cells (p53-) were arrested preferentially 
in G2/M phase, which is indicative of non-functional p53 
(Figure lb, bottom panels). Importantly, cell cycle profiles 
of both cell lines were similar in the absence of DNA 
damage (Figure lb, upper panels). Thus, collectively. 



our data confirmed that U2-0S cells, in which p53 was stably 
knocked down by specific shRNA behaved similarly to 
authentic p53-null cells. 

To uncover specific microRNAs regulated by p53 during cell 
cycle arrest in response to doxorubicin, a global analysis of 
microRNA expression patterns in U2-0S scr and KDp53 cells 
was performed (data not shown). Among several microRNAs 
that were regulated by p53 (including miR-34a), we focused 
on two highly expressed microRNAs: miR-26a-1 and miR-1 6-2. 
As a validation of our microarray results Q-RT-PCR of 
miR-16 and miR-26a showed a statistically significant 
increase of their expression levels (1.4 and 1.6 fold, 
respectively) on 24 h of doxorubicin treatment versus control 
cells (Figures 1c and d, respectively). Notably, the level of 
miR-1 6 and miR-26a expression in U2-0S KDp53 cells versus 
the wild-type U2-0S cells was lower even in the absence of 
DNA damage, suggesting that p53 regulates these genes at 
basal level. Importantly, DNA damage enhanced miR-16 and 
miR-26a expression in U2-0S KDp53 cells at a lower rate 
compared with U2-0S scr cells with wild-type p53 (2.2 and 4 
fold for miR-16 and miR-26a, respectively). As a positive 
control for p53 function, we analysed transcription of the p21 
gene, which is a known p53 target in these cells (Figure 1e). 
Collectively, our data strongly support the results of micro- 
array experiment, which identified miR-26a and miR-16 as 
p53- and doxorubicin-dependent microRNAs. 

p53 controls transcription of miR-1 6-2 and miR-26a-1 in 
luciferase assay. It should be noted that the mature 
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Figure 1 Expression of miR-16 and miR-26a is induced in the presence of p53 by doxorubicin, (a) Comparison of the p53 levels in U2-0S cells stably expressing 
scrambled shRNA versus p53-specific shRNA (U2-0S KDp53) on treatment with doxorubicin. U2-0S cells with wild type and knocked-down expression of p53 were 
continuously treated with 0.5 of doxorubicin for 0, 6, and 12 h before harvesting. Cell lysates were then prepared and analysed by western blotting using p53-speicific 
antibody. Anti-tubulin serum was used as loading control, (b) Cell cycle distribution of U2-0S cells with wild type and knocked-down expression of p53 on doxorubicin 
treatment. Both cell types were treated with 0.5 of doxorubicin for 12 h, fixed with formaldehyde and stained with propidium iodine for subsequent cell cycle analysis. 
The same cells were treated with doxorubicin for 24 h and total RNA was extracted, converted to cDNA and analysed by real-time PCR for expression levels of miR-16 
(c); miR-26a (d), and p21, as a positive control (e) 
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products of miR-16 and 26a each comprise RNAs tran- 
scribed from two different genes. Specifically, miR-16 is 
encoded by the miR-16-1 and miR-16-2 genes located on the 
chromosomes 4 and 13, respectively. Analogously, miR-26a 
is the product of expression of two genes, miR-26a-1 and 
miR-26a-2, located on the chromosome 3 and 14, respectively. 
Thus, to elucidate which specific genes contribute to the 
upregulation of miR-16 and miR-26a by p53, we tested 
promoters of the respective genes in the luciferase reporter 
assay (Figures 2a and b). We found that miR-16-1 and miR- 
26a-2 promoters were not dependent on the presence of p53 
(Figure 2b and data not shown, respectively). Both in the 
absence or presence of doxorubicin treatment the luciferase 
activity of the miR-26a full-length promoter was higher than 
that of the deletion mutant lacking several putative p53- 
binding sites (region from -1600 to -1400) (Figure 2a). 
Importantly, in U2-0S KDp53 cells the transcriptional activity 
of both miR-26a promoter constructs was significantly 
attenuated compared with the cells with wild-type p53 even 
after doxorubicin treatment (Figure 2a, right panel). 

Similar results were obtained for the miR-16-2 promoter 
(Figure 2b). Three fragments encompassing 3200 bp of 
the upstream sequence were tested individually for their 
response to doxorubicin and p53. Fragments 1 and 2 of the 
miR-16-2 promoter exhibited the highest activity in the 
p53-positive cells both in the absence and in the presence 
of doxorubicin (Figure 2b, left panel and right panel. 



respectively). Importantly, ablation of p53 in U2-0S KDp53 
cells correlated with a threefold decrease of luciferase activity 
for both miR-16-2 fragment 1- and 2-containing constructs 
(Figure 2b, right panel). The miR-16-1 promoter-containing 
construct (a kind gift of G Calin) displayed the activity only 
marginally higher than the background level, indicating that 
this promoter did not depend on the presence of p53. Taken 
together, these data confirm the results of microarray 
experiments suggesting that both miR-26a-1 and miR-16-2 
are upregulated by p53 in response to doxorubicin. 

p53 directly binds to the regulatory regions of miR-16-2 
and miR-26a-1 genes. To determine whether p53 was 
directly bound to the promoters of miR-16-2 and miR-26a-1 
genes, we performed ChIP assay using p53-specific anti- 
body (Figure 3).To independently validate the results of 
luciferase assay, we analysed p53 binding in the region from 
-2000 to +3000 of the miR-26a-1 gene using a series of 
primers (Figure 3a). As shown in Figure 3b, maximal p53 
binding after DNA damage was detected in the amplicon 
-1500. Also, p53 binding was detected in the amplicon 
-2000. Importantly, much lower binding of p53 was 
observed at the amplicons -500 and +3000, indicating 
that the maximal p53 binding occurs in the region from 
-2000 to -1500. p53 binding was significantly lower in 
U2-0S KDp53 cells, thus confirming the specificity of p53 
binding in our assay. 
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Figure 2 p53 and doxorubicin treatment activate promoters of miR-1 6-2 and miR-26a genes in luciferase assay, (a) Shown is the scheme of the promoter and its deletion 
mutant for the miR-26a-1 gene situated in front of the luciferase gene. U2-0S cells stably expressing scrambled shRNA (U2-0S scr) versus p53-specific shRNA (U2-0S 
KDp53) were treated or not treated with doxorubicin for 14 h before harvesting. Cell extracts were then prepared, normalized for transfection efficiency by measuring beta- 
galactosidase activity and assayed for luciferase activity. Empty luciferase vector with minimal HSP70 promoter was used as negative control, (b) Shown are three fragments 
of the miR-16-2 promoter, which were assayed for their ability to drive luciferase activity in response to doxorubicin. Cells were treated and assayed as described above. 
A fragment of the miR-16-1 promoter driving transcription of the luciferase gene that does not respond to p53 was used as a negative control 
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Figure 3 Analysis of p53 binding to tine regulatory regions of miR-16-2 and miR-26a genes by ChIP assay, (a) Shown is the scheme of the miR-26a-1 regulatory 
region. Putative p53-binding sites are depicted. Positions of nucleotides are given in respect to the transcription start site, (b) A summary of results of ChIP assay for p53 
biniding on U2-0S scr (WT) and U2-0S shRNA-p53 (KD) cells. Cells were non-treated or treated for 24 h with doxorubicin and analysed for p53 binding to different 
amplicons of the sequence of the miR-26a-1 gene, (c) ChIP assay of Drosha and (d) p68 binding to the region - 1500, where maximal p53 binding was detected, 
(e) Shown is the scheme of the SMC4 gene with intronic regulatory region for miR-16-2. (f) ChIP assay of p53 binding in the promoter of SMC4 gene and (G) in the 
upstream sequence of miR-16-2 



As p53 was shown to interact with the Drosha micro- 
processing compiex,^^'^"^ we hypothesized that p53 can 
recruit components of the Drosha complex to the target 
promoters thus enhancing the processing of specific miRs. To 
test this hypothesis, we performed ChIP assay using 
antibodies against Drosha and p68, which is an ATPase 
enzyme often associated with Drosha.^"^ The maximal binding 
of Drosha was detected at the replicon -1500 of the 
miR-26a-1 promoter in the p53-specific manner (Figure 3c). 
Unexpectedly, the recruitment of p68 did not depend on the 
presence of p53, but did depend on the treatment with 
doxorubicin (Figure 3d). Thus, we found that Drosha, but not 
p68, is recruited by p53 to the promoter of miR-26a-1 on DNA 
damage. 

We also tested the ability of p53 to bind the regulatory 
regions of miR-16-2 (Figure 3e). Our luciferase reporter 
experiments identified the region -2000 to + 1 as the most 
responsive to p53. Thus, to analyse p53 binding we used an 
amplicon -700 bp, which is located between the promoter 
fragments 1 and 2. It should be noted that miR-16-2 is coded 
by an intronic gene, located within the SMC4 (structural 
maintenance of chromosomes) gene (Figure 3e). The 
promoter region of SMC4 has been shown to contain 
the binding sites for E2F1 to control miR-16-2/miR-15b 
expression. Since E2F1 binding sites often overlap with 
the p53 ones, we decided to test whether p53 also binds 
this region upon DNA damage ( - 4 kb). As shown in Figures 
3f and g, p53 bound to both regions ( - 4 kb and - 700 bp) in 
the time- and doxorubicin-dependent manner reaching its 
maximum at 12h. 



p53 regulates the expression of miR-16-2 and miR-26a-1 
by different meclianisms. Previously published studies 
suggested that p53 has role in regulation of microRNA 
expression both on the transcriptional and post-transcrip- 
tional levels. Having this in mind, we examined the 
expression of miR-16 and miR-26a both at the level of 
transcription (primary miRs) and after their processing 
(mature miRs) (Figures 4a and b, respectively). Analysis of 
the primary miRs (pri-miRs) levels suggests that only the 
miR-26a-1 gene, but not miR-16-1 or miR-16-2, was 
regulated transcriptionally by p53 in response to doxorubicin 
(Figure 4a). On the contrary, both mature miR-26a and 
miR-16 showed their dependency on p53 and doxorubicin 
(Figure 4b). Collectively, these results suggest that p53 
controls the expression of miR-26a and miR-16 differently, 
that is, miR-26a transcriptionally and miR-16 post-transcrip- 
tionally, on the level of its maturation. 

miR-16 and miR-26a target checlcpoint Icinases Clield 
and Weel . Next, we wished to identify the targets of miR-16 
and miR-26a. miR-16 has been reported to target several cell 
cycle regulators, including cyclin E, Chkl, cyclin D1, and 
Cdk6.^^'^^ To identify the targets for miR-26a, we compared 
the lists of the first 50 putative targets predicted by PicTar 
and MiRanda algorithms (http://pictar.mdc-berlin.de/ and 
http://www.microrna.org, respectively). Among a few com- 
mon targets found in these lists, we chose Weel for future 
experiments, due to its pharmacological relevance. Importantly, 
Weel was predicted to be also the target for miR-16 
(Figure 5a). To test whether miR-26a and miR-16 indeed 
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Figure 4 p53 differentially affects expression of primary and mature miR-16-2 
and miR-26a. Effect of p53 activation by doxorubicin on the transcription of primary 
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target Chk1 and Wee1, we employed a luciferase reporter 
assay in which the activity of the luciferase gene fused to 
3'UTR of either Wee1 or Chk1 was assayed in the presence 
of miR-16 or miR-26, or both together (Figure 5b). Transfec- 
tion of miR-16 decreased the Luc 3-UTR Weel construct 
activity compared with control, suggesting that miR-16 
directly targets Weel. Unexpectedly, despite the bioinfor- 
matics prediction, miR-26a did not affect the activity of 
Luc 3-UTR Weel (Figure 5b). However, when both miRs were 
co-transfected together, a synergistic effect was observed, 
resulting in a more than threefold decrease of the luciferase 
activity. Similarly, activity of the Luc 3-UTR Chekl construct 
was also attenuated fivefold versus the control activity when 
both miRs were co-transfected (Figure 5b). 

Next, we decided to confirm the previously published 
targets for miR-16 by transfecting synthetic pre-curser of 
miR-16 or its inhibitor, antago-miR-16, into p53-deficient 
HI 299 cells and probe the levels of cyclin E and Chkl 
expression by western blotting (Figure 5c). As evident from 
the figure, 15nM of miR-16 pre-cursor completely eliminated 
cyclin E and severely attenuated the level of Chkl expression. 
This effect was specific, because the expression levels of both 
cyclin E and Chkl were restored when miR-16 pre-cursor was 
co-transfected with antago-miR-16 inhibitor (Figure 5c). 

To further validate our results obtained in the luciferase test, 
we examined how miR-16 and miR-26a pre-cursors affect 
Weel and Chkl on the protein level (Figures 5d and e). 



On transfection with miR-26a pre-cursor, the level of Weel did 
not change, which is consistent with the results of the 
luciferase assay (Figure 5d, upper panel). Unexpectedly, 
despite the fact that miR-26a has no predicted binding sites in 
the 3'-UTR of Chkl, it affected the Chkl protein levels 
(Figure 5d, middle panel). Consistent with the results of 
the luciferase assay, ectopic expression of miR-16 attenuated 
the expression of both Weel and Chkl (twofold). Notably, the 
levels of Weel and Chkl expression were synergistically 
downregulated (threefold and fivefold, respectively) by 
co-transfection of the miR-16 and miR-26a pre-cursors 
(Figure 5d, middle panel). 

We hypothesized that as DNA damage activates expres- 
sion of miR-16 and miR-26 mediated by p53, then the 
expression of Chkl and Weel should be negatively regulated 
by DNA damage in a p53-dependent manner. To validate this 
hypothesis, isogenic HCT116 cell lines (p53+ and p53") 
were untreated or treated with doxorubicin and analysed for 
the presence of Chkl and Weel proteins (Figure 5e). 
In accord with our predictions, the levels of both Chkl and 
Weel proteins were diminished on activation of p53 in 
HCT116 p53+ cells, but not in the p53~ ones (Figure 5e). 
Taken together, these results suggest that expression of the 
Weel and Chkl proteins may be controlled by p53-dependent 
miR-16 and -26a. 

Overexpression of miR-16 and miR-26a on DNA damage 
augments G1/S cell cycle arrest and apoptosis. As both 
Chkl and Weel kinases participate in regulation of the cell 
cycle progression through G2/M, we sought to ascertain their 
effect on cell cycle on DNA damage by miR-dependent 
ablation of these proteins (Figure 6). U2-0S cells transfected 
with scrambled miR, miR-16, or miR-26a were non-treated or 
treated with doxorubicin. In the absence of doxorubicin, all 
three cell lines showed similar cell cycle distribution patterns 
(Figure 6a, upper panel). However, doxorubicin treatment 
affected these cell lines differently: control cells were 
arrested in G2/M phase (60%) with only 20% arrested in 
G1 and 5% undergone apoptosis (Figure 6a, lower panel, 
and Figure 6b), whereas cells expressing miR-16 or miR-26a 
exhibited a more pronounced cell cycle arrest in G1 phase 
(30% and 35%, respectively) and an increased apoptosis 
(20% and 15%, respectively). 

Augmented expression of Chekl and Wee1 is associated 
with poor survival prognosis. To assess the biological 
significance of our in vitro findings on Weel and Chkl 
inhibition by p53-dependent miR-16 and miR-26a, we 
determined the role of these checkpoint kinases in survival 
of breast cancer patients using bioinformatics (PPISURV,^^ 
bioprofiling.de^^) (Figure 7, upper panels). Both kinases had 
negative impact on the survival of breast cancer patients. 
Next, we wanted to determine whether the in vitro effects of 
miR-16 and miR-26a correlated with survival rates of patients 
diagnosed with breast or prostate cancer using a publically 
available bioinformatics tool MiruMir^° (Figure 7). Importantly, 
high levels of expression for miR-16 and miR-26a correlated 
with better survival rates for both breast and prostate cancer 
patients. These results suggest that miR-16 and miR-26a are 
important prognostic markers and that their targets, including 
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Figure 5 Chk1 and Wee1 are the targets of p53-dependent miR-16 and miR-26a. (a) Effect of miR-16 overexpression alone or together with its antago-mlR (inhibitor) 
on expression levels of Chkl and cyclin E was assessed by western blotting using specific antibodies. Beta-tubulin was used as a loading control, (b) Effect of individual 
or simultaneous overexpression of miR-16 and miR-26a on expression levels of Weel and Chkl in HI 299 cells. Expression levels of Weel and Chkl in the cells transfected 
with control oligo were arbitrarily set as 1 . Intensities of bands corresponding to the Weel and Chkl proteins were calculated using ImageQuant software, (c) To determine the 
effect of p53 activation on cellular levels of the Chkl and Weel proteins isogenic HCT1 16 cell lines with different p53 status were treated with doxorubicin for 16 h followed by 
western blotting against the respective proteins. Note, that anti-p53 antibody gave a non-specific band in samples prepared from p53-deficient cells, (d) Shown are the regions 
of homology between Weel and Chkl 3-UTRs and miR-1 6 and miR-26a sequences, (e) A schematic presentation of putative miR-1 6- and miR-26a-binding sites in the 3-UTR 
regions of Chkl and Weel . Below are shown the results of luciferase assay in which the activities of 3-UTRs of Chkl and Weel fused to the luciferase gene constructs were 
measured in the presence of individually or simultaneously expressed miR-16 and/or miR-26a. Binding of these mlRs to 3-UTRs of the respective genes should result in 
repression of the protein synthesis and hence decrease luciferase activity. Fold repression was calculated against the value of luciferase activity of the cells transfected with 
respective vectors in the absence of mlRs, which was set as 1 



Wee1 and Chk1 , negatively contribute to the survival rates of 
cancer patients. 

We also investigated whether these two microRNAs were 
co-regulated in breast cancer samples and whether this 
phenomenon was dependent on p53. To this end, we 
computed individual Pearson correlations for miR-16 and 
miR-26a expression profiles with all other measured 410 
microRNAs from 51 different breast cancer cell lines 
described in Riaz et al^^ The value of correlation between 
miR-16 and miR-26a was 0.57 and has rank 7 for both 
microRNAs among 410 microRNAs tested. This provides an 
estimation of significance of the correlation: the probability to 
select randomly a partner for either miR-16 or miR-26a with 
correlation higher than 0.57 is calculated as 6/41 0~ 0.014. 
This calculated P-value confirms that the correlation between 
miR-16 and miR-26 is significant. 

Moreover, when considering only the TP53-positive cell 
lines, Pearson correlation for both miR-26a and miR-16 
expression profiles increased to 0.71, which statistically 
supports our hypothesis that TP53 is closely involved in 
regulation of both microRNAs. 



Discussion 

The biological response to DNA damage requires the 
activation of the cell death pathway by p53 or its family 
members. ^^"^"^ p53 is a powerful transcription factor^^ that 
drives a large number of promoters, including novel 
pathways such as the connection between IL-7Ra and 
telomere erosion^^ or the silencing of repeats and non-coding 
RNA,^^ and the metabolism of the cell via the pathways of 
mevalonate^^ or serine.'^^ As a consequence of transactiva- 
tion of so many different promoters described so far, p53 is 
able to control different biochemical pathways affecting the 
regulation of life'^'' and death"^^ of the cell, or its senescence"^^ 
and hypoxia."^"^ 

Because of such versatility, p53 is considered as a critical 
tumour suppressor and thus, perhaps not surprisingly, found 
mutated or deleted in more than half of all human tumours. 

Our microarray data showed that among several 
microRNAs regulated by p53 in response to genotoxic stress, 
miR-16 and miR-26a belonged to the most abundantly 
expressed ones. Importantly, both luciferase assay and ChIP 
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Figure 6 Cell cycle analysis of U2-0S cells transfected with scramble, miR-16, 
and miR-26a oligonucleotides, (a) Cell cycle profiles of U2-0S cells transfected with 
indicated miR oligos in the absence or presence of 0.5 /aU doxorubicin. Following 
transfection, the cells were treated with doxorubicin for 14 h before fixing with 
methanol and staining with propidium iodine, (b) Graphic representation of the 
results shown in a 



analysis showed that p53 bound the promoters of miR-16 and 
miR-26a genes (Figures 2 and 3). However, despite the 
activation of miR-16-2 promoter observed in the luciferase 
assay, we failed to detect any appreciable transcriptional 
induction of the endogenous miR-16-2 gene. This result 
argues that the artificial reporter systems do not always 
faithfully reproduce regulatory events on the endogenous 
promoters. Importantly, we also found that p53 regulated 
miR-16-2 not transcriptionally, but post-transcriptionally, 
which is consistent with previously published data.^^ 

The p53-regulated miR-16 and miR-26a are known to 
have anti-proliferative and pro-apoptotic effects in cancer 
cells. Significantly, the miR-26a-1 -encoding locus at 
3p21 .3 is contained within a region that is frequently deleted in 
small cell lung carcinomas, renal cell carcinomas, and breast 
carcinomas, suggesting that this micro-RNA is a tumour 
suppressor."^^ Likewise, the miR-1 6-1 -containing locus is 
either deleted or inactivated in chronic lymphocytic leukaemia, 
prostate cancer, and non-small lung carcinoma.^^"^^ Among 
the known targets of miR-16 are the cell cycle-regulating and 
anti-apoptotic proteins (cyclin D1 and D3, cyclin E, cdk6, and 
bcl-2, respectively).^^'^^ 

Our bioinformatics analysis predicted that in addition to the 
targets mentioned above, miR-16 and miR-26a may also 
target cell cycle kinases, Chkl and Weel. In line with our 
findings is a recently published report showing that the 
miR-1 5/1 6/1 95/424/497 family sensitized cisplatin-resistant 
cells to apoptosis by targeting the WEE1 and CI-IK1 gene 



expression. Our data suggest that miR-26a by itself did 
not affect the translation levels of 3-UTR'Wee1-luc, but 
synergistically inhibited its expression when co-transfected 
with miR-16. This finding was confirmed by western blot 
analysis (Figure 5b). It can be hypothesized that when acting 
together, these miRs de-stabilize the 3'-UTR of Weel more 
efficiently to prevent the loading of translation factors. 
Although no binding sites for miR-26a was predicted in the 
3'-UTR of the CIHEKI gene, we detected a 50% decrease 
of the expression level of Chkl on ectopic expression of 
miR-26a, as judged by western blotting (Figure 5b). One 
possibility is that miR-26a targets the coding region of Chkl, 
which is unusual for miRs. Alternatively, miR-26a may affect 
the expression of Chkl indirectly, by targeting its positive 
regulator(s). 

The significance of our findings is underscored by the 
fact that both Chkl and Weel are important pharmacological 
targets. As a consequence of their inactivation, cancer 
cells undergo unscheduled mitosis and halt in G1 or get 
eliminated via apoptosis. Accordingly, genetic ablation or 
pharmacological inhibition of Chkl and/or Weel results in 
sensitization of cancer cells to genotoxic therapy.^^ Moreover, 
several clinical studies suggest that the pharmacological 
inhibition of these kinases is a new promising therapy to treat 
various tumours. 

The biological importance of our results is supported by the 
bioinformatics analysis across several expression platforms, 
whereby the expression of both miR-16 and miR-26a 
correlates with a better survival prognosis for breast cancer 
and prostate cancer patients (Figure 7). 

In addition, our data provide a new insight to a long-standing 
observation that p53 represses Weel and Chkl on genotoxic 
stress.^"^'^^ Thus, according to our hypothesis, upregulation of 
miR-16 and miR-26a by p53 during genotoxic stress 
attenuates the expression of Chkl and Weel kinases, 
resulting in the loss of G2/M checkpoint and unscheduled 
mitosis. This causes an accumulation of cells in G1 phase and 
their apoptosis (Figure 8). 

Materials and Methods 

Cells. The human lung adenocarcinoma cell line H1299, human osteosarcoma 
cell line U2-0S, and human embryonic kidney cell line HEK-293 were purchased 
from the ATCC (IVIanassas, VA, USA). Human colon adenocarcinoma cell lines 
HCT1 1 6 (p53 + ) and p53( - ) were kindly provided by B Vogelstein. U2-0S cells 
with stable expression of shRNA against p53 were obtained by lentiviral 
transfection with either empty or bearing shRNA against p53 pLKO vectors as 
described previously.^^ These cells were maintained in DMEM medium (Gibco, 
Carlsbad, CA, USA) or DMEM medium (Gibco) supplemented with 10% fetal 
bovine serum (Gibco). Cells were grown at 37 °C in a humidified atmosphere with 
5% CO2. 

miR-16 and miR-26a overexpression or knockdown. miR-16 and 
miR-26a overexpression was achieved by transfecting cells with pre-miR-16 
(a synthetic RNA oligonucleotide duplex mimicking miR-16 precursor), whereas 
miR-16 knockdown was achieved by transfecting cells with anti-miR-16 
(a chemically modified single-stranded antisense oligonucleotide designed to 
specifically target against mature miR-16). Scrambled negative control RNA 
(pre-miR-control and anti-mlR-control) served as negative control. Synthetic RNA 
molecules, including pre-miR-16, anti-miR-16 and scrambled negative control 
RNA, were purchased from Dharmacon (Thermo Scientific). Cells were seeded on 
60-mm dishes and were transfected the following day using Lipofectamine 2000 
(Invitrogen, Carlsbad, CA, USA) according to the manufacturer's instructions. 
For each well, equal doses (200pmol) of pre-miR-control, pre-miR-16, anti-mi R-control 
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Figure 7 A bioinformatics analysis of the effect of Chkl , Weel , miR-16 and miR-26a expression on survival rates for both breast and prostate cancer patients. The two 
upper panels show correlations between the levels of expression of Chkl and Weel and survival rates of breast cancer patients. The two middle and two bottom panels show 
correlations between miR-16 and miR-26a levels of expression and the survival rates of breast and prostate cancer patients, respectively 
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Figure 8 A schematic model of how p53 via activation of miR-26a and miR-16 
promotes G1/S arrest and apoptosis in response to genotoxic drugs. On the one 
hand, miR-16 and 26a downregulate G2/M checkpoint kinases, Chkl and Weel, 
resulting in slippage of damaged cancer cells from the mitotic checkpoint arrest and 
their elimination via apoptosis. On the other hand, miR-16 represses expression of 
several cell cycle promoting genes, including cyclin E, cyclin D, cdk4 and cdk6 
resulting in their arrest in G1/S, which subsequently elicits apoptosis 



or anti-miR-16 were added. Similar procedures were done for miR-26a. Cells were 
harvested 24 h after transfection. 

RNA isolation and relative quantification RT-PCR. Total RNA was 
extracted from the cultured cells using TRIzol Reagent (Invitrogen) according to 
the manufacturer's instructions. Assays to quantify mature miR-16 and miR-26a 



were carried out using microRNA kit (Applied Biosystems, Foster City, CA, USA) 
according to the manufacturer's instructions. miRNA expression in cells was 
normalized to that of the U6 snRNA. The amount of miR-16 relative to the internal 
control U6 was calculated with the equation 2~^^^\ in which AACt = 

(Cj miR — Cj ueltarget — (Cj miR-16 — Oj uelcontrol- 

For relative quantification RT-PCR analysis of p21 and j6-actin mRNA, 1 fig of 
total RNA was reverse-transcribed to cDNA with oligod(T) using Ready-to-Go kit 
(GE Healthcare, Amersham, UK). Then real-time PCR was performed on 
RotorGene 6000 PCR machine (Qiagen, Hilden, Germany) using SYBR green 
dye (Invitrogen). The 20-|il PCR reaction included 1 ^1 RT product, 1 x QuantiTect 
SYBR green PCR Master Mix, and 0.5 ^M each sense and antisense primers. The 
reactions were incubated in a 96-well plate at 95 °C for 5 min, followed by 40 cycles 
of 95 °C for 30 s, 60 °C for 30 s, and 72 °C for 30 s. All reactions were run in 
triplicate. After the reactions, the Cj values were determined using fixed threshold 
settings. The relative amount of p21 mRNA was normalized to p-acWn mRNA. 
Expression of p21 was analysed by RT-PCR using the following primers: 5'-G 
ACACCACTGGAGGGTGACT, 3'-CTCTTGGAGAAGATCAGCCG. Primers to 
beta-actin were purchased from (PrimerDesign, Southampton, UK). 

Microarray procedures. RNA was extracted by using TRIzol (Invitrogen), 
according to the manufacturer's protocol. RNA was labelled and hybridized on 
microRNA microarray chips as previously described.^^ Briefly, 5 ^g of RNA from 
each sample was biotin labelled during reverse transcription using random 
hexamers. Raw data were normalized and analysed by GeneSpring GX software 
version 7.3 (Agilent Technologies, Santa Clara, CA, USA) as described.^^ 

miR-16 and miR-26a target prediction. The miRNA target 
prediction and analysis was performed with the algorithms from TargetScan 
(http://www.targetscan.org/) PicTar (http://pictar.mdc-berlin.de/) and miRanda 
(http://www.microrna.org/). 
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Chip assay. Chromatin immunoprecipitation (ChIP) assays were performed as 
described previously.^^ Briefly, 3x10^ cells per sample were cross-linked with 1% 
formaldehyde for 15 min. Cross-linking was neutralized with 0.125 M glycine, and cells 
were scraped in PBS. Chromatin was sonicated using the Diagenode (Liege, 
Belgium) Bioruptor for 15 min with 30-s pulse/pause cycles in polycarbonate tubes 
on ice to shear chromatin to 300- to 600-bp fragments. Unsheared debris was 
spun down, and then the chromatin was incubated overnight with the appropriate 
antisera, concurrent with the blocking of protein G-Sepharose beads using 2.5% 
BSA. Immune complexes were then precipitated using 'blocked' protein G beads 
for 4h at 4°C, washed three times, and then eluted. Immunoprecipitated DNA 
was purified and quantitative PCR was performed with 1 ^1 of DNA to assess 
p53, Drosha and p68 binding levels. The following antibodies were used: 
anti-p53 (Ab-6, Millipore, Billerica, MA, USA), anti-Drosha (Abeam, Cambridge, 
UK; ab37165) and anti-p68 (Santa Cruz Biotechnology, Santa Cruz, CA, USA). 
The following primers were used to analyse miR-26a: +500 bp: 5'-C 
CTTCTTTCAGTCACCCATGA, 3'-GCTGGCCGAGGAATTCTATAC; +1500 bp, 
5'-CTTCAACCCCAGACCTTCCT; 3'-TCCAGGCCCAGCATATTG; +2000 bp, 5'-C 
TGCCATTGCAGGTGTATGT, 3'-GCCAAGCTGGCAGTAAAGAG; -3000 bp, 5'-TG 
GATCCCAAGTGAAGAACC; 3-GAAAAACCATTCCCCAGCTT. 

CHIP primers for mlR 15b/16-2: SMC4 promoter/lst exon: 5'-GCAGGAGCGA 
CAATAAGATGG; 3'-CGCTCCTACCGGTGTTTCG; miR-16-2 (-700) (700 bp 
up from the miR 15b/16-2 gene, p53BS) 5'-TTCCATGCTGTTTTGGTCCC, 
3'-AGGCTTGTTTGGAGGTGGGT. 

Western blotting. p53, Chkl, Cyclin E and Weel protein levels were 
quantified by western blot analysis of whole-cell extracts using antibodies against 
the corresponding proteins. These samples were normalized by blotting with an 
antibody against a-tubulin. The following antibodies were used: anti-p53 (Ab-6, 
Millipore), anti-Cyclin E (HE12; Santa Cruz Biotechnology), anti-Weel (Cell 
Signalling, Danvers, MA, USA), Chkl (Cell Signalling), and anti-a-tubulin (B-7, 
Santa Cruz Biotechnology). 

Luciferase assay of 3'-UTR constructs. The entire 3'-UTRs of 
human Weel and Chekl were amplified from human genomic DNA using PCR. 
The PCR products were inserted as Xbal fragments into the p-MIR-report 
plasmid (Ambion, Life Technologies, Carlsbad, CA, USA). The following primers 
were used to amplify Weel 3'UTR: 5'-CCAAGTCTAGAGAACCGCTCTGTC 
AGCCTTA and 3'-AAAGATCTAGAAGAAAACAGGCATCACAAGGA. Chkl 3'UTR 
was amplified using the following primers: 5'-ACCGATCTAGAAGCCAGAAG 
ATTTGGCTTCC and 3'-TTCGGTCTAGAGGCTTCGCTTCACAGACTGA. 

Efficiency of the cloning was confirmed by sequencing. For luciferase reporter 
assays, cells were cultured in six-well plates, and each well was transfected with 
0.5 lag of firefly luciferase reporter plasmid, 0.5 /ug of j6-galactosidase expression 
plasmid (Ambion), and equal amounts of scrambled negative control RNA, 
pre-miR-26a, pre-miR-16, or anti-miR-16 using Lipofectamine 2000 (Invitrogen). 
The j6-galactosidase plasmid was used as a transfection control. At 24 h post 
transfection, cells were assayed using luciferase assay kits (Promega, Madison, Wl, 
USA). The data depicted are representative of three independent experiments 
performed on different days. 

Cell cycle assay. Cells were harvested, washed once with PBS, and fixed in 
70% ethanol overnight. Staining for DNA content was performed with 50 mg ml" 
propidium iodide and 1 mg ml" ^ RNase A for 30 min. Analysis was performed on 
a FACS Canto II flow cytometer (BD Biosciences, Franklin Lakes, NJ, USA) with 
Cell Quest Pro software. Cell cycle modelling was performed with Modfit 3.0 
software (Verity Software House, Topsham, ME, USA). 

Statistical analysis. All presented images of western blotting and cell cycle 
assay are representative of at least three independent experiments. Relative 
quantification RT-PCR, luciferase reporter, and cell viability assays were 
performed in triplicate, and each experiment was repeated three to five times. 
The data shown are the mean ± S.D. of at least three independent experiments. 
Statistical significance was considered at P<0.05 using the Student's f-test. 

mlRs expression analysis in respect to survival outcome. Effects 
of miR-16 and miR-26a expression on survival outcome of cancer patients were 
analysed using a bioinformatics tool, MIRUMIR. This algorithm exploits rank 
information from the publically available clinical miR expression data sets 
annotated with patient survival information. For each available data set, samples 



are grouped with respect to expression rank of the user-specified miR. The 'low 
expression' and 'high expression' groups are those where expression rank of miR is 
less or more than average expression rank across the data set, respectively. This 
separation of patients into 'low' and 'high' groups along with survival information is 
next used to find any statistical differences in survival outcome. 
The R statistical package is used to perform survival analyses and to draw KM plots. 
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